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0\ ; ABSTRACT 

I discuss the decay modes of the neutral and charged Higgs bosons in the minimal 
$H , supersymmetric extension of the Standard Model. Special emphasis will be put 

on the the QCD corrections to the hadronic decay modes, the below threshold 
-three body- decays and the decays into supersymmetric particles, charginos, 
neutralinos and sfermions. A Fortran code calculating the various Higgs decay 
branching ratios is then briefly presented. 



1. Introduction 
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In the Minimal Supersymmetric extension of the Standard Model 1 (MSSM), the 
Higgs sector 2 is extended to comprise three neutral h/H (CP=+), A (CP=-) and a 
pair of charged scalar particles H ^ . The Higgs sector is highly constrained since there 
are only two free parameters at tree-level: a Higgs mass parameter [generally Ma] and 
the ratio of the vacuum expectation values of the two doublet fields responsible for 
the symmetry breaking, tg/3 [which in Grand Unified Supersymmetric models with 



d.: 

QJ , b-r Yukawa coupling unification 3 is forced to be either small, tg/3 ~ 1.5, or large, 

tg/3 ~ 50]. After the inclusion of the large radiative corrections 4 , while the lightest 
Higgs boson h is predicted to be lighter than M h < 130 GeV, the H, A and H ± states 
^ ■ are expected to have masses of the order of a few hundred GeV. 

The decay pattern of the MSSM Higgs bosons is determined to a large extent by 
their couplings to fermions and gauge bosons, which in general depend strongly on 
tg/3 and the mixing angle a in the CP-even sector. The pseudoscalar and charged 
Higgs boson couplings to down (up) type fermions are (inversely) proportional to tg/3; 
the pseudoscalar A has no tree level couplings to gauge bosons. For the CP-even 
Higgs bosons, the couplings to down (up) type fermions are enhanced (suppressed) 
compared to the SM Higgs couplings [tg/3 > 1]; the couplings to gauge bosons are 
suppressed by sin / cos(/3 — a) factors [see Table 1.] 

For large values of tg/3 the pattern is simple, a result of the strong enhancement of 
the Higgs couplings to down-type fermions. The neutral Higgs bosons will decay into 
bb (~ 90%) and r + r~ (~ 10%) pairs, and into rv T pairs below and tb pairs above 
the top-bottom threshold. For the CP-even Higgs bosons h and H, only when Mh 

*Talk given at the Ringberg Workshop the Higgs Puzzle, Ringberg Castle, Tegernsee, Germany, 
December 8-13 1996; to appear in the proceedings. 
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approaches its maximal value is this simple rule modified: in this decoupling limit, 
the h boson is SM-like and decays into charm and gluons with a rate similar to the 
one for t + t~ [~ 5%] and in the high mass range, Mh ~ 130 GeV, into W pairs with 
one of the W bosons being virtual; the H boson will mainly decay into hh and AA 
final states. 





9&UU 


5W<2 




h 


cos a/ sin f3 


-> 1 


— sin a/ cos /5 


-> 1 


sin(/3 — a) — > 1 


H 


sin a/ sin (3 — > 


l/tg/3 


cos a/ cos (5 — > 


tg/3 


cos(/3 — a) — > 


A 


l/tg/3 




tg/3 








Table 1: Higgs couplings to fermions and gauge bosons normalized to the SM Higgs cou- 
plings, and their limit for Ma 3> Mz- 



For small values of tg/3 ~ 1 the decay pattern of the heavy neutral Higgs bosons 
is much more complicated. The b decays are in general not dominant any more; 
instead, cascade decays to pairs of light Higgs bosons and mixed pairs of Higgs and 
gauge bosons are important and decays to WW/ZZ pairs will play a role. For very 
large masses, they decay almost exclusively to top quark pairs. The decay pattern of 
the charged Higgs bosons for small tg/3 is similar to that at large tg/3 except in the 
intermediate mass range where cascade decays to Wh are dominant. 

When the decays into supersymmetric particles are kinematically allowed, as it 
should be the case at least for the heavy CP-even, CP-odd and charged Higgs bosons, 
the pattern becomes even more complicated since the decay channels into charginos, 
neutralinos and squarks play will play a non-negligible role. 

In the following, I will discuss three topics related to the decay modes of the Higgs 
particles in the MSSM: (a) the QCD corrections to the hadronic decay modes 5 , (b) 
the below threshold three-body decays 6 and (c) the decays into SUSY particles 7 of 
the heavy H, A and H ± bosons, including the QCD corrections to the squark decay 
modes 8 . I will then briefly introduce a Fortran code 9 which calculates the various 
decay branching ratios. For more details and for a complete list of references, see the 
original papers Refs. [5-9]. 



2. Hadronic Decay Modes: QCD Corrections 5 



The particle width for decays to massless b, c quarks directly coupled to the Higgs 
particle is given, up to 0(a 2 s ) QCD corrections [the effect of the electroweak radiative 
corrections in the branching ratios is negligible], by the well-known expression 



r[* - QQ] = 3 -^gi QQ rnUMz) 



1 + 5.67— + (35.94 - 1.36JV»-£ 

TT TT Z 



(1) 



in the MS renormalization scheme; the running quark mass and the QCD coupling 
are defined at the scale of the Higgs mass, absorbing this way any large logarithms. 
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The quark masses can be neglected in general except for top quark decays where this 
approximation holds only sufficiently far above threshold. Since the relation between 
the charm pole mass M c and the MS mass evaluated at the pole mass m c (M c ) is 
badly convergent, one can adopt the running quark masses TFTqIMq) [which have been 
extracted directly from QCD sum rules evaluated in a consistent 0(a s ) expansion] as 
starting points. The evolution from Mq to a scale fx ~ M$ is given by: 

m Q (fi) = m Q (M Q )c[a s (fi)/7r}/c[a s (M Q )/7T} 



c( x ) = (25/6x) 12/25 [1 + 1.014a: + 1.39x 2 } for M c < fx < M b 
c{x) = (23/6x) 12/23 [1 + 1.175a; + 1.50a; 2 ] for M b < fx 



(2) 



Typical values of the running b, c masses at the scale fx = 100 GeV, characteristic 
for M$, are displayed in Table 2, with the evolution calculated for a s (M z ) = 0.118 ± 
0.006; Mq 2 are the quark pole masses. 





a s {M z ) 


Mq(Mq) Mq = Mq Z 


m Q {fi -- 100 GeV) 


b 


0.112 
0.118 
0.124 


(4.26 ± 0.02) GeV (4.62 ± 0.02) GeV 
(4.23 ± 0.02) GeV (4.62 ± 0.02) GeV 
(4.19 ± 0.02) GeV (4.62 ± 0.02) GeV 


(3.04 ±0.02) GeV 
(2.92 ±0.02) GeV 
(2.80 ±0.02) GeV 


c 


0.112 
0.118 
0.124 


(1.25 ± 0.03) GeV (1.42 ± 0.03) GeV 
(1.23 ± 0.03) GeV (1.42 ± 0.03) GeV 
(1.19 ± 0.03) GeV (1.42 ± 0.03) GeV 


(0.69 ±0.02) GeV 
(0.62 ±0.02) GeV 
(0.53 ±0.02) GeV 



Table 2: The running b and c quark masses in the MS scheme at a scale fi = 100 GeV. 

The decay of the Higgs bosons to gluons is, to a good approximation mediated by 
heavy top quark loops; the partial decay width, including QCD radiative corrections 
which are built up by the exchange of virtual gluons and the splitting of a gluon into 
two gluons or into Np massless quark-antiquark pairs, is given by [fx ~ M$] 



r N »[<t>^99 + -] 
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(3) 



16^3 

with <p = h,H and a s = a^ F (fi 2 ). The radiative corrections are very large, nearly 
doubling the partial width. The final states $ — > bbg and ccg are also generated 
through processes in which the b, c quarks are coupled to the Higgs boson directly. 
Gluon splitting g —> bb in $ — > gg increases the inclusive decay probabilities r($ — > 
bb ± . . .) etc. Since b quarks, and eventually c quarks, can in principle be tagged 
experimentally, it is physically meaningful to consider the particle width of Higgs 
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decays to gluon and light u, d, s quark final jets separately. The contribution of 6, c 
quark final states to the coefficient in front of a s in eq. (3) is: 



7 1 M| M | 
-3 + 3 [l ° g M! + hg M! ] 

Instead of naively subtracting this contribution, it may be noticed that the mass 
logarithms can be absorbed by changing the number of active flavors from Np = 5 
to N F = 3 in the QCD coupling a^ Np h The subtracted parts may be added to the 
partial decay widths into c and b quarks. 

The numerical analysis of the branching ratios for the lightest CP-even Higgs 
decays in the decoupling limit where h is SM-like, with the quark masses and QCD 
couplings given above and a top mass M t = (176 ±11) GeV, is shown in Fig. 1. To 
estimate systematic uncertainties, the variation of the c mass has been stretched over 
2a and the uncertainty of the b mass to 0.05 GeV. However, the dominant error in 
the predictions is due to the uncertainty in a s and the errors in the prediction for the 
charm and gluon branching ratios are very large. Nevertheless, the expected hierarchy 
of the Higgs decay modes is clearly visible despite these uncertainties. Similar results 
hold for the heavy CP-even and CP-odd Higgs decays. 
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Fig. 1: Branching ratios of the h boson in the decoupling limit, including the uncertainties 
from the quark masses and the QCD coupling a s [shaded bands]. 



4 



3. Three Body decay modes 6 

Besides these two-body decays, below-threshold modes can play an important 
role. It is well-known that SM Higgs decays into real and virtual Z pairs are quite 
substantial: the suppression by the off-shell propagator and the additional Zff cou- 
pling is at least partly compensated by the large Higgs coupling to the Z bosons. 
For the same reason, three-body decays of MSSM Higgs particles mediated by gauge 
bosons, heavy Higgs bosons and top quarks, are of physical interest. Important three- 
body decays for the H, A and H ± bosons are [V = W, Z\: 

H -> VV* -> VfP , AZ* -> Aff , H ± W T * -> if*//' , it* -> ibW + (4) 
A -> hZ* -> /i// , tt* -> tW+ (5) 
fi* -> -> hff , AH/* -> A//' , It* -> &&lf (6) 

For the lightest Higgs boson /i, the only releveant below threshold decay mode is 
h -> W*W* for M fe ~ 130 GeV. In this case, both the W's have to be taken off-shell. 
The branching ratios for h, H, A and H ± decays are shown in Fig. 2 for tg/3 = 1.5, in 
the case where the mixing in the stop sector is neglected. 

For the heavy Higgs boson if, the decay H — > hh is the dominant channel, 
superseded by ti decays above the threshold [for the latter, the inclusion of the three- 
body modes provides a smooth transition from below to above threshold]. This rule 
is only broken for Higgs masses of about 140 GeV where an accidentally small value 
of the Xnhh coupling allows the bb and WW* decay modes to become dominant. 
Important channels in general, below the ti threshold, are decays to pairs of gauge 
bosons and bb decays. In a restricted range of M H , below-threshold AZ* and H ± W Zf * 
also play a non-negligible role. In the case of the pseudoscalar A, the dominant modes 
are the A — > bb and A — > ti decays below the hZ and ti thresholds respectively; in the 
intermediate mass region, Ma = 200 to 300 GeV, the decay A — > hZ* [which reaches 
~ 1% already at Ma = 130 GeV] dominates. The gluonic decays are significant 
around the ti threshold. For the charged Higgs boson, the inclusion of the three-body 
decay modes will reduce the branching ratio for the tu channel quite significantly. 
Indeed, this decay does not overwhelm all the other modes since the three-body decay 
channels H + — > hW* as well as H + — ► AW* in the low mass range and H + — > bt* in 
the intermediate mass range have appreciable branching ratios. 

The total widths of the Higgs bosons are in general considerably smaller than for 
the SM Higgs due to the absence or the suppression of the decays to W/Z bosons 
which grow as M%. The dominant decays are built-up by top quarks so that the 
widths rise only linearly with M$. However, for large tg/3 values, the decay widths 
scale in general like tg 2 /9 and can become experimentally significant, for tg[3 > (9(30) 
and for large M$. 
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Fig. 2: Branching ratios for the CP-even, the CP-odd and the charged MSSM Higgs bosons, 
including the three-body decays, for tg/3 = 1.5 and no stop mixing. 
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4. SUSY Decay modes 7 



In the previous discussion, we have assumed that decay channels into neutralinos, 
charginos and sfermions are shut. However, these channels could play a significant 
role, since some of these particles can have masses in the 0(100 GeV) range or 
less. To discuss these decays, we will restrict ourselves to the MSSM constrained 
by minimal Supergravity, in which the SUSY sector is described in terms of five 
universal parameters at the GUT scale: the common scalar mass m , the common 
gaugino mass M1/2, the trilinear coupling A, the bilinear coupling B and the higgsino 
mass fj,. These parameters evolve according to the RGEs, forming the supersymmetric 
particle spectrum at low energy. 

The requirement of radiative electroweak symmetry breaking further constrains 
the SUSY spectrum, since the minimization of the one-loop Higgs potential specifies 
the parameter \i [to within a sign] and also B. The unification of the b and r Yukawa 
couplings gives another constraint: in the A t fixed-point region, the value of tg/3 is 
fixed by the top quark mass through: m t ~ (200 GeV) sin/3, leading to tg/3 ~ 1.75. 
There also exists a high-tg/3 [A& and A T fixed-point] region for which tg(3 ~ 50-60. 
If one also notes that moderate values of the trilinear coupling A have little effect 
on the resulting spectrum, then the whole SUSY spectrum will be a function of tg/3 
which we take to be tg/3 = 1.75 and 50, the sign of /x, m which in practice we replace 
with Ma taking the two illustrative values Ma = 300 and 600 GeV, and the common 
gaugino mass M ± / 2 that we will freely vary. 

The decay widths of the heavy CP-even, the CP-odd and the charged Higgs 
bosons, H, A and H^, into pairs of neutralinos and charginos [dashed lines], squarks 
[long-dashed lines] and sleptons [dot-dashed lines], as well as the total [solid lines] 
and non-SUSY [dotted-lines] decay widths, are shown in Fig. 3 for tg/3 = 1.75, /i > 
and two values of M A = 300 [left curves] and 600 GeV [right curves]. 

For Ma = 300 GeV, i.e. below the ti threshold, the widths of the H decays 
into inos and sfermions are much larger than the non-SUSY decays. In particular, 
squark [in fact t and b only] decays are almost two-orders of magnitude larger when 
kinematically allowed. The situation changes dramatically for larger Ma when the 
ti channel opens up: only the decays into t pairs when allowed are competitive with 
the dominant H — > ti channel. Nevertheless, the decays into inos are still substantial 
having BRs at the level of 20%; the decays into sleptons never exceed a few percent. 

In the case of the pseudoscalar A, because of CP-invariance and the fact that 
sfermion mixing is small except in the stop sector, only the decays into inos and 
A — > i\i 2 decays are allowed. For these channels, the situation is quite similar to 
the case of H: below the ti threshold the decay width into ino pairs is much larger 
than the non-SUSY decay widths [here t 2 is too heavy for the A — > i\t 2 decay to be 
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allowed], but above 2m t only the A — > tit 2 channel competes with the tt decays. 

For the charged Higgs boson H^, only the decay H + — > t±bi [when kinematically 
allowed] competes with the dominant H + — > tb mode, yet the x + X° decays have a 
branching ratio of a few ten percent; the decays into sleptons are at most of 0(%). 

In the case where fj, < 0, the situation is quite similar as above. For large tg/3 
values, tg/3 ~ 50, all gauginos and sfermions are very heavy and therefore kinemati- 
cally inaccessible, except for the lightest neutralino and the r slepton. Moreover, the 
bb/rr and tb/ris [for the neutral and charged Higgs bosons respectively] are enhanced 
so strongly, that they leave no chance for the SUSY decay modes to be significant. 
Therefore, for large tg/3, the simple pattern of bb/rr and tb decays for heavy neutral 
and charged Higgs bosons still holds true even when the SUSY decays are allowed. 
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Fig. 3: Decay widths for the SUSY decay modes of the heavy CP-even, CP-odd and charged 
Higgs bosons, for tg/3 = 1.75. The total and the non-SUSY widths are also shown. 
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Since the decays into stop and sbottom squarks can be dominant when kinemati- 
cally allowed, QCD corrections must be incorporated in order to have full control on 
the decay widths and to make a reliable comparison with the standard (non-SUSY) 
decay channels. The QCD corrections to the decays of the heavy CP-even, CP-odd 
and charged MSSM Higgs bosons into stop and sbottom quarks have been recently 
calculated 8 . These corrections are found to be rather large, enhancing or supressing 
the widths by amounts up to 50% and in some case even more. The QCD correc- 
tions depend strongly on the gluino mass; however for very heavy gluinos, they are 
only logarithmically dependent on nig. Contrary to the case of Higgs boson decays 
into light quark pairs, these large corrections cannot be absorbed into running squark 
masses since the latter are expected to be of the same order as as the Higgs masses. 

5. The program HDECAY 9 

Finally, let me shortly describe the fortran code HDECAY, which calculates the 
various decay widths and the branching ratios of Higgs bosons in the SM and the 
MSSM and which includes: 

(a) All decay channels that are kinematically allowed and which have branching 
ratios larger than 10~ 4 , y compris the loop mediated, the three body decay modes 
and in the MSSM the cascade and the supersymmetric decay channels. 

(b) All relevant two-loop QCD corrections to the decays into quark pairs and to 
the quark loop mediated decays into gluons are incorporated in the most complete 
form; the small leading electroweak corrections are also included. 

(c) Double off-shell decays of the CP-even Higgs bosons into massive gauge bosons 
which then decay into four massless fermions, and all all important below-threshold 
three-body decays discussed previously. 

(d) In the MSSM, the complete radiative corrections in the effective potential 
approach with full mixing in the stop/sbottom sectors; it uses the renormalisation 
group improved values of the Higgs masses and couplings and the relevant leading 
next-to-leading-order corrections are also implemented. 

(e) In the MSSM, all the decays into SUSY particles (neutralinos, charginos, 
sleptons and squarks including mixing in the stop, sbottom and stau sectors) when 
they are kinematically allowed. The SUSY particles are also included in the loop 
mediated 77 and gg decay channels. 

The basic input parameters, fermion and gauge boson masses and total widths, 
coupling constants and in the MSSM, soft-SUSY breaking parameters can be chosen 
from an input file. In this file several flags allow to switch on/off or change some 
options [e.g. chose a particular Higgs boson, include/exclude the multi-body or 
SUSY decays, or include/exclude specific higher-order QCD corrections]. The results 
for the many decay branching ratios and the total decay widths are written to several 
output files with headers indicating the processes and giving the input parameters. 
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